Model-independent tracking of criticality signals in nuclear multifragmentation data 
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We look for signals of criticality in multifragment production in heavy-ion collisions using model- 
independent universal fluctuations theory. The phenomenon is studied as a function of system size, 
bombarding energy, and impact parameter in a wide range of INDRA data. For very central collisions 
{b/bmax < 0.1) we flnd evidence that the largest fragment in each event, Zmax, plays the role of 
an order parameter, deflning two different regimes at low and high incident energy, respectively, 
according to the scaling properties of its fluctuations. Data for a wide range of system masses and 
incident energies collapse on to an approximately universal scaling function in each regime for the 
most central collisions. The forms of the scaling functions for the two regimes are established, and 
their dependence on the total mass and the bombarding energy is mapped out. Data suggest that 
these regimes are linked to the disappearance of heavy residues in central collisions. 



I. INTRODUCTION 

It has long been hoped that the study of nuclear multi- 
fragmentation reactions as observed in intermediate en- 
ergy heavy-ion collisions Q, |^ |E| can give valu- 
able information on the nuclear matter phase diagram 
and equation of state 0, S S| ■ The principal guide- wire 
for this research has been the search for signs of some- 
thing analogous to a liquid-gas phase transition in data 
on intermediate mass fragment (IMF) production, ever 
since the observation of power laws in fragment mass- 
yield distributions 0, ^^li reminiscent of the condensa- 
tion of liquid drops in a critical vapour . This effort 
was encouraged by the failure of statistical models of hot 
nuclear decay Il2. fl^ to explain the observed frag- 
ment yields [2, llal unless they suppose the break up of 
the system at low densities |T3, Ha] where the phase 
transition may occur 

Signals of the phase transition in experimental data 
on multifragmentation may be revealed by anomalously 



large fluctuations of fragment observables 0, US l2ll |. 
The main obstacles to such endeavours insofar as exper- 
imental data are concerned are the huge statistical fluc- 
tuations inherent to small systems such as atomic nuclei, 
and the still-open question of the mechanism(s) of frag- 
ment production in heavy ion collisions at intermediate 
enerries^^2i BI2I iMlpliliip^illlllip, 
IH^ l3(iL I37I l38L ■ Indeed in all experimental studies 
of the question, at the same time as one is searching for 
evidence of the phase transition of a piece of hot nuclear 
matter one is (implicitly or not) obliged to eluci- 
date the manner in which such an excited sj;stem may be 
formed in the course of certain reactions [3J|- The solid- 
ity of any experimental evidence for a link between mul- 
tifragmentation and the nuclear matter phase diagram 
will be undermined by any remaining ambiguity about 
the dominant mechanism of fragment production in the 
selected sample of experimental events, and it would be 
preferable to be able to address the question of "phases" 
in nuclear multifragmentation independently of reaction 
mechanism. 



* deceased 



It is for this reason that the theory of the universal 
character of order parameter fluctuations in flnite sys- 
tems proposed by Botet and Ploszajczak |4(l lilL 142. li^ 
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provides an attractive opportunity to address the ques- 
tion of phase transitions induced by heavy-ion coUisions 
in the least ambiguous way possible. According to their 
work, it is possible to obtain pertinent information on the 
relationship between the formation of clusters in a system 
and the phase transition (s) of said system without need- 
ing to characterise the state of the system under study to 
the extent of, for example, supposing it to be in thermo- 
dynamical equilibrium at the time of cluster formation. 
One needs only to study the properties of the clusters: 
specifically, all pertinent information can be obtained 
from a sufHciently precise measurement of the event-by- 
event distributions of cluster multiplicity and the size of 
the largest cluster produced. The attractiveness of such 
an approach in the domain of nuclear fragmentation re- 
actions is obvious, for the reasons given above: although 
the mechanism of fragment production remains an open 
question, for several years now large solid-angle multide- 
tector arrays such as INDRA have provided high-quality 
data on the multiplicity and the size (charge) of the frag- 
ments produced in such reactions. 

A first application of the universal fiuctuations the- 
ory to INDRA data for central collisions of Xe+Sn from 
25 to 50 A MeV was pubHshed in 0|- In this paper we 
will present the results of the same analysis applied to a 
wide range of systems measured with INDRA. As in 
we observe a signature compatible with the existence of 
different regimes at different bombarding energies. This 
behaviour disappears for less central collisions. By apply- 
ing the analysis to data for colliding systems with total 
mass number between 73 and 394 we will show the mass 
dependence of the energy of transition between the two 
regimes, as well as presenting a systematic study of the 
universal scaling functions observed in the low- and high- 
energy regimes. 



II. UNIVERSAL FLUCTUATIONS OF THE 
ORDER PARAMETER IN FINITE SYSTEMS 

Universal scaling laws of fluctuations (the A-scaling 
laws) have been derived for equilibrium systems, and 
shown to apply also in certain out-of-equilibrium situ- 
ations. In a system in which the second-order critical 
behaviour can be identifled, the relation between order 
parameter, criticality and scaling law of fluctuations has 
been established and the relation between the scaling 
function and the critical exponents has been found. De- 
tails can be found in (i^ . 

Experimental observables that may be related to a crit- 
ical order parameter can be identifled through their A- 
scaling behaviour. The A-scaling is observed when two 
or more probability distributions P/v [m] of the observable 
m for a system of 'size' N collapse onto a single scaHng 
curve <i>(2;(A)) independent of system size when plotted 
in terms of the scaling variables: 



<m>^ P^H = Hha)) = ( "^^^y ) (1) 

where < m > is the mean value of the distribution Pn [m] 
and ^<A<1. <m> plays the role of a scale param- 
eter and can replace TV as a measure of the size of the 
system. A less strong (necessary but not suflicient) condi- 
tion for A-scaling is that the variance of the distribution 
should scale with its mean value as 



a' ^< m (2) 

so that in a log- log plot of cr^ versus < m >^ data should 
fall on a straight line of slope A. 

The scaHng law Equation with A = ^ is associated 
with low temperature ("ordered") systems, or with ob- 
servables which are not related to an order parameter. 
ScaHng with A = 1 is seen at high temperature ("dis- 
ordered" system) and also for critical systems. For m 
to be an order parameter it must exhibit a correspond- 
ing change of A-scaling regime when some suitable con- 
trol parameter (e.g. available energy, temperature, bond- 
breaking probability, etc.) is varied. 

Here it is worth saying a word about the comparison of 
experimental distributions using EqQJand the determina- 
tion of A-scaling regimes for data. As the transformation 
from experimental observable to the scaling variable 
is a linear transformation, the form of the <i>(2:(A)) dis- 
tribution is the same as that of PAr[m]. However, the 
presence of the exponent A in the scale factor < to >^ 
means that two identical distributions will appear differ- 
ent (their widths wiH be different) if they are scaled using 
a value of A which is not the one relating the mean values 
and variances of the two distributions via Eq[2 For ex- 
ample, consider two Gaussian distributions whose widths 
and mean values are related via EqlJlwith A = 1. In this 
case the use of EqQlwith, for example, A = ^ would not 
lead to a universal scaling function: the widths of the 
two scaled distributions would be different, and, because 
of the normalisation of <i>(2:(A)), the height of the distri- 
butions would differ also. It should also be noted that 
the relationship between the mean and the variance of 
a set of distributions, i.e. the value of A if Eq[2l holds, 
is quite independent of the form of the distributions: a 
given value of A does not imply a certain type of distri- 
bution, and vice versa. 

More detailed information on the state of the system 
and the fragmentation process may be found in the form 
of the scaling functions <i>(2:(A))) Eq. ^ For systems 
far from a critical point, the central limit theorem tells 
us that for an observable m being the sum of uncorre- 
lated random variables, one should observe asymptoti- 
cally a Gaussian distribution for fluctuations of to about 
its mean value. If on the other hand ra is an extremal 
value such as the largest among a set of uncorrelated 
random variables, then asymptotically its distribution 
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should be that of Gumbel's first asymptote (i^. For 
critical systems such general results do not exist, as, due 
to the presence of correlations at all length-scales, the 
order parameter distribution must depend on the precise 
details of the interaction in this case. An asymptotic 
form of the large-m scaHng function tail was derived in 
[i^ . where close to a critical point the order parameter 
distribution was expected to fall off like exp —m'^ with 
3 < i> < 6, meaning that large deviations of the order 
parameter from its mean value are strongly suppressed. 

There are two generic families of fragment produc- 
tion scenarios for which the second-order phase transition 
has been identified, with two different order parameters. 
These are 

• the fragment multiplicity for fragmentation scenar- 
ios such as the fragmentation-inactivation binary 
(FIB) model 

• the size of the largest cluster or fragment in the so- 
called aggregation scenarios such as percolation or 
Fisher droplet models [4Si |. 

Therefore the elimination of one of these two easily mea- 
surable experimental quantities as not having the order 
parameter A-scaling behaviour described above should 
give important information on the fragment production 
process, by allowing to exclude one of the two scenarios. 

III. PRESENTATION OF DATA 
A. Experimental details 

In order to study as exhaustively as possible the ques- 
tion of the existence of an order parameter or other 
"phase-like" behaviour in heavy-ion collisions, we have 
profited from the wide range of very high quality data 
which has been obtained with the INDRA Air array 
El 19] at the GANIL (Caen) and GSI (Darmstadt) 
accelerator facilities. This charged product detector cov- 
ers about 90% of the 47r solid angle. The total number 
of detection cells is 336 arranged according to 17 rings 
centred on the beam axis. The first ring (2° — 3°) was 
composed of 12 fast NE102/NE115 phoswich detectors 
during the experiments at the GANIL facility. For the 
GSI experiments these were replaced by 12 telescopes 
composed of a 300 /im thick silicon detector and a 14 cm 
thick CsI(Tl) scintillator. Rings 2 to 9 cover the angular 
range from 3° to 45° and are made of three detector lay- 
ers: a low pressure gas-ionization chamber (5cm of CsFg 
at 20 to 50 mbar, depending on the experiment and the 
polar angle), a 300 ^m thick silicon detector and a 14 to 
10 cm thick CsI(Tl) scintillator. The remaining 8 rings 
cover the angular range from 45° to 176° and have two de- 
tection layers: ionization chamber and 7.6 to 5 cm thick 
CsI(Tl) scintillator. Fragments with Z up to the charge 
of the projectile are identified with unit resolution in the 
forward region, when they are stopped in the scintillator 



detectors. Beyond 45°, the charge resolution is one unit 
up to ^ = 16, and a few charges for larger Z. Over the 
whole angular range, a very good isotope identification 
is obtained for Z = 1 to Z = 3, except for particles with 
low energies where ambiguities are unresolved. 

The energy calibration of the CsI(Tl) scintillators was 
obtained for light charged particles (LCP) by means of 
the elastic and inelastic scattering of secondary LCP 
beams (p,d,t,'^He,^He) produced by the fragmentation of 
a 95 A MeV beam in a thick C target. These par- 
ticles were then momentum selected by the "alpha mag- 
netic spectrometer" of GANIL and scattered in a C or 
Ta target installed in the INDRA reaction chamber. A 
typical energy resolution was about 4%. Typical identi- 
fication thresholds are a few 100 keV for light particles, 
0.7^ MeV for Z = 3, and lAA MeV for Z = 35. A com- 
plete technical description of INDRA its calibration and 
its electronics can be found in (ij lia l49llH(l IHlI IHIIH^ . 



B. Overview of data 
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Figure 1: (Color online) Experimental correlations between 
the size of the largest detected fragment Zmax and the to- 
tal transverse energy of light charged particles divided by the 
available cm. energy, for Xe+Sn collisions. A minimum of 
4 well-identified charged particles was required in the off-line 
analysis. Logarithmic contour levels are shown corresponding 
to the number of events (darker tones indicate larger num- 
bers). 

The data presented here cover a wide range of quasi- 
symmetric systems studied with the INDRA array, with 
different total masses and bombarding energies. They 
are: 
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Ar+KCl 32, 40, 52 and 74A MeV; 

Ni+Ni 32, 40, 52, 63, 74, 82 and 90A MeV; 

Xe+Sn 25, 32, 39, 45, 50, 65*, 80* and 100*A MeV; 

Au+Au 40*, 60* and 80*A MeV. 

The systems marked with an asterisk were measured dur- 
ing the experimental campaign at the GSI facihty. We 
have concentrated on symmetric colliding systems in or- 
der to benefit from the maximum overall efficiency of the 
INDRA array in this case. Before presenting the analy- 
sis of this data set in terms of universal ffuctuations we 
will give an overview of the chief characteristics of these 
reactions. 

Figure n] shows experimental data for the Xe+Sn sys- 
tem. The contour plots show the number of events 
measured corresponding to each value of the size of the 
largest detected fragment, Z^ax, and of the fraction of 
the available energy converted into transverse energy of 
light charged particles, Eti2 / Egyau . This latter quantity 
has been shown (23,l^l5J;l^ to be principally related 
to the geometry of heavy-ion collisions in this energy do- 
main, and is particularly well-suited to sorting events 
measured with the INDRA detector with little bias, be- 
cause the efficiency of the array for light charged parti- 
cle detection is « 90% whatever the centrality /reaction 
mechanism. The data shown were recorded with an on- 
line trigger requiring that at least four detectors fired in 
coincidence, while in the off-line analysis we required at 
least four correctly identified charged products in each 
considered event. Very similar plots were obtained for 
all other data studied in this paper, as for example the 
data for the Ni+Ni system shown in Figure [21 They give 
an overview of the evolution of reaction mechanisms with 
beam energy and impact parameter. 

For the least violent collisions (small Eti2/ Eavau) two 
strong contributions to the total cross-section are ob- 
served, depending on whether the projectile-Hke fragment 



was detected (Z„ 



-'pro] , 



or not [Zmax < 10). For 



most of these collisions the target-like residue is too slow- 
moving to be detected and/or correctly identified by the 
INDRA array. With increasing beam energy, the contri- 
bution from projectile-like fragments appears to decrease 
in importance. This is due to the increasingly forward- 
focused kinematics of the reactions and the smaller graz- 
ing angle at higher energy, which mean that projectile- 
like fragments from peripheral collisions are more likely 
to be defiected to angles too small for them to be detected 
in the first ring of the INDRA array fsee lTTmi . 

For more central collisions (larger Eti2 / Eavaii) the out- 
come of the reaction depends on bombarding energy and 
the mass of the coUiding nuclei. For the Xe+Sn sys- 
tem at the lowest incident energies (Figures QJi-b) , the 
cross-section for all but the most peripheral collisions is 
dominated by events for which the heaviest detected frag- 
ment has around one half of the charge of the projectile 
nucleus, and there is a discontinuous evolution in cross 




Figure 2: (Color online) Experimental correlations between 
the size of the largest detected fragment Zmax and the total 
transverse energy of light charged particles divided by avail- 
able cm. energy, for Ni+Ni collisions. A minimum of 4 well- 
identified charged particles was required in the off-line anal- 
ysis. Logarithmic contour levels are shown corresponding to 
the number of events (darker tones indicate larger numbers). 



section between the most peripheral collisions and these 
events. This is due to the onset of fission of projectile-like 
nuclei above a certain threshold excitation energy which 
is attained in mid-peripheral and more central reactions 
[55| . At higher bombarding energies (Figures Qt-d) a 
more continuous evolution of projectile-like residue size 
with the coUision violence is seen. For a fighter system 
such as Ni+Ni (Figures |2K-b) a continuous ridge goes 
from the most peripheral toward more central collisions 
showing that the size of the projectile-like residues de- 
crease continuously with decreasing impact parameter, 
as fission is not a predominant decay channel even at 
high excitation energies for such light nuclei. The cor- 
relations between Zmax and Etu/ Eavaii for the Ar+KCl 
system (not shown) are very similar to those for Ni+Ni, 
whereas for the Au+Au system they resemble those for 
Xe+Sn at 50 and lOOA MeV except that for the most 
peripheral collisions a clear contribution from the fission 
of gold quasiprojectiles is seen. 

Finally let us remark the similarity between Figures 
and [2t, concerning the production of heavy residues 
at low incident energy. Both figures show that there are 
two distinct contributions to the total yield for residues 
with Z close to that of the projectile. The first was dis- 
cussed in the previous paragraph and is due to projectile- 
like fragments produced in peripheral collisions (small 
Eti2/ Eavaii values). This is the dominant contribution. 
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Nevertheless there is a second contribution which is as- 
sociated with a broad distribution of large residue sizes 
Zproj/'^ < Zrnax ^ ^proj as wcll as an equally broad dis- 
tribution of Eti2/ Eavaii values corresponding to mid- to 
central collisions. The appearance of these experimental 
correlations suggests that the mean value of the charge of 
the heaviest detected fragment, < Z^ax >, first decreases 
with increasing collision violence and then increases for 
the most "central" collisions at the lowest beam energy 



one considers all (Figure El left panel) or only well- 
measured (right panel - ratio of total detected charge, 
Zdet, to total charge of projectile and target required 
to be at least 80%) events. The requirement of well- 
measured events reveals the monotonous decrease of the 
mean charge of projectile-like fragments in peripheral col- 
lisions {Etu/Eavaii < 0.1) by excluding from the average 
those events where neither target-like nor projectile-like 
fragments were detected (compare Figure^. 
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Figure 3: (Color online) The mean charge of the largest frag- 



ment detected in each event, Zrr 



as a function of the total 



transverse energy of light charged particles (normalised to the 
available centre of mass energy) for Xe+Sn collisions between 
25 and 50A MeV. Vertical bars show the estimated statisti- 
cal error on the mean. It should be noted that events with 
Etu/ Eavaii > 0.2 correspond to only 1% of the total mea- 
sured cross-section. 



for these two systems. 

According to Figures and El such reactions make a 
relatively important contribution to the production of 
heavy residues in central collisions of Ni+Ni and Xe+Sn 
at 32 and 25A MeV respectively. When the incident en- 
ergy increases (40^ MeV for Ni+Ni, FigureEt, and 
MeV for Xe+Sn, FigureQlD) the size and yield of the asso- 
ciated residues decrease making this contribution harder 
to distinguish. From Figures Et-d and Figures Qt-d it 
appears that the cross-sections associated with these re- 
actions become very small above ~ 4QA MeV, for the 
Ni+Ni system, or above ^ 32A MeV, for the Xe+Sn sys- 
tem (for the Ar+KCl system, not shown here, this con- 
tribution is discernible up to h2A MeV). However, more 
sensitive analyses (e.g. see below for Xe+Sn) can reveal 
the survival of such reactions at higher incident energies. 

Figure El shows the evolution of < Zmax > with colli- 
sion centrality for the Xe+Sn system between 25 and 50^ 
MeV. It can be seen that the mean charge of the heaviest 
fragment produced in each event increases with increas- 
ing centrality for 25, 32, and, very sHghtly, S9A MeV. It 
should be noted that this tendency is observed whether 




Figure 4: (Color online) Collisions from the most peripheral 
with detection of a quasi-projectile residue (furthest right- 
hand side of the figure) to the most central (indicated by a 
ring) are classed into bins corresponding to 1% of the total 
measured cross section. 



Another way of excluding poorly-measured events is 
to keep only those for which most of the momentum of 
the incident beam particles is reconstructed from the de- 
tected nuclei. This gate on data has the effect of retaining 
all events for which at least a good reconstruction of the 
quasi-projectile residue and products was obtained, for 
the most peripheral collisions, whilst in the most central 
collisions we achieve an almost complete reconstruction 
of the event. In this way we can follow, in Figure ^ 
the correlated evolution of the variance and the mean of 
the Zmax distribution as a function of collision violence: 
from peripheral collisions leading to a slightly excited 
quasi-projectile (large < Zmax > and small variance) to 
the most central coUisions (indicated by a ring around 
the last data point for each beam energy). Each point 
in this figure corresponds to Eti2 cuts defined by slicing 
the minimum-bias distribution into 100 bins, each con- 
taining an equal number of events. Thus each point rep- 
resents 1% of the measured cross-section, although due 
to our requirement of well-measured events the most pe- 
ripheral (least well measured) events are absent. It allows 
to see how in the most central coUisions two different be- 
haviours are observed depending on the incident beam 
energy: below 39A MeV in the 3-5% most central col- 
lisions the mean charge of the largest fragment actually 
begins to increase again, whilst the correlation between 
this mean value and the variance of the distribution is not 
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the same as for the majority of other impact parameters. 

This figure suggests that at beam energies ^39^ 
MeV the origin of the heaviest fragment is not the same 
in central collisions as for the rest of the reactions. The 
increase of < Z„mx > for the most central collisions at 
these energies suggests a contribution from (incomplete) 
fusion reactions where the heaviest fragment detected is 
an evaporation residue. The disappearance of this phe- 
nomenon for >39A MeV can be interpreted as signalling 
the disappearance of fusion-evaporation residues, and 
marking the onset either of fusion-multifragmentation 
or of transparency in central collisions (incomplete 
stopping) 



transparen 
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IV. 



ANALYSIS OF DATA IN TERMS OF 
UNIVERSAL FLUCTUATIONS 

A. Selection of central collisions 



We begin our analysis by extending the results of [43 |. 
for central collisions of Xe+Sn, to a wider range of sys- 
tems measured with the INDRA multidetector array. As 
was discussed in Sec. HI the analysis of universal fluc- 
tuations in multifragmentation data does not require to 
know (or to suppose) how fragments are formed. However 
the comparison of events with very different collision ge- 
ometries is unlikely to give meaningful results. Therefore 
we will limit our study to very central collisions for which 
the geometrical overlap between projectile and target is 
as close as possible to total. The need to select equiva- 
lent classes of events for a large range of system masses 
and beam energies with as little dependence as possible 
on detector efficiency led us to use the total transverse 
energy of light charged particles, Eti2- 

The sorting variable EtA 2 has been studied specifically 
for the INDRA detector |HE El ■ It was shown that, for 
a given projectile-target system, the minimum-bias in- 
clusive distributions of Eti2 for different beam energies 
scale as a function of the available centre of mass energy, 
consistent with Etu being mainly sensitive to the geom- 
etry of the collisions, i.e. that it is a good indicator of 
collision centrality. Moreover, as it was pointed out in 
Sec. 



InTBl the efficiency of the INDRA detector for fight 
charged particles is almost independent of the type of re- 
action under study, allowing to use this variable to sort 
all data considered in this paper in the same way. 

In addition, in this way we avoid any direct link be- 
tween the studied observables (size of the largest frag- 
ment in each event and fragment multiplicity) and the 
variable used for the selection of events. We are there- 
fore in the best possible situation to avoid distortion of 
the data due to autocorrelations with the experimental 
filter [H. 

The data analysed in the following with the largest 
values of Efu correspond to 1% of the total number of 
events recorded during the experiment with a minimum- 
bias condition (6 < 0-lbmax in the geometrical approx- 



imation of where the maximum impact parameter 
bmax is smaller than the sum of the radii of projectile and 
target due to the experimental trigger condition). We 
were able to check whether the largest detected charged 
fragment of each event is really the largest, using the 
total detected charge of the event, Ztot- 



B. Establishing an "order parameter" for nuclear 
multifragmentation with no model-dependent 
hypotheses 
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Figure 5: (Color online) Data for the total multiplicity of 
charged products with Z > 1, Mtot, for well- measured central 

collisions of Xe + Sn (6 < O.lbmax, Zdet > 0.8{Zproj + Ztarg))- 

(a) Log-log plot of the variance versus the squared mean value 
of the distribution of Mtot for each bombarding energy. The 
straight line represents Eq. |2|with A = i. (b) Experimental 
distributions of Mtot expressed in the variables of the second 
scaling law, i.e. Equation U with A — ^. The dashed curve is 
a Gaussian fit to all the data. 

Generic models of cluster production may be classed 
into two types. The first, the class of "fragmentation" 
scenarios in which a system is broken up by a series of 
binary splittings or some other physical process, has for 
its order parameter the number or multiplicity of clus- 
ters. For such models, it is the multiplicity which ex- 
hibits different A-scafing regimes if the system has differ- 
ent phases. The second class of models, in which clusters 
are built up by "aggregation" of smaller constituents, has 
for its order parameter the size of the largest cluster. As 
INDRA does not measure the mass but only the atomic 
number of fragments, we will assume in the following that 
the largest fragment of each event corresponds to Z^ax, 
the fragment with the largest atomic number. 

Figure El shows the data for the total multipficity of 
charged products with Z > 1, Mtot, for central collisions 
of Xe+Sn from 25 to lOOA MeV. In FigureEt we plot the 
natural logarithm of the variance of the measured Mtot 
distributions as a function of the natural logarithm of the 
square of the mean value. The estimated statistical errors 
of these quantities are smaller than the symbols used. If 
fluctuations obey the universal scaling law Eq. [2l then 
the data must fall on a straight line of slope A in this 
plot. It can be seen in FigureEt that this is true to a fair 
approximation, and that the variance of the multiplicity 
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distribution grows with increasing bombarding energy as 
< Mtot >. 

Figure Eh shows that the multipUcity distributions for 
different bombarding energies collapse to a unique distri- 
bution (scaling function) when expressed in terms of the 
second scaling law (Equation^with A = i), as suggested 
by Figure Et- It should be noted that examination of an 
observable's scaling properties in this way is far more 
constraining than that of Fig. I^i, when the statistics 
of the data samples allow it. The observed scaHng func- 
tion is very well approximated by a Gaussian distribution 
(dashed curve on the figure). The scaling properties of 
total multiplicity fiuctuations are therefore the same for 
all bombarding energies in the range 25-100^ MeV. This 
is true not only for the Xe+Sn data but also for all the 
data we have studied in this paper: see for example the 
data for ^®Ni+^®Ni collisions in Figure[2K- Therefore the 
total multiplicity of charged particles for central colli- 
sions in this energy range does not show any evidence of 
'anomalous' or 'critical' behaviour. 
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Figure 6: (Color online) As Figure|^but now for the multiplic- 
ity of fragments with Z > 3, Mz>3. (b) Data for bombarding 
energies 32-65^ MeV expressed in the variables of the second 
scaling law, i.e. Equation U with A = i. 

The multiplicity Mtot is dominated by the multiplicity 
of light charged particles, Mlcp, which is typically 3-4 
times greater than the multiplicity of IMF {Z > 3). How- 
ever, similar conclusions as for Mtot can be drawn for the 
IMF multiplicity, Mz>3, for which the Xe+Sn data are 
presented in FigureEl the Ni+Ni data in FigureEjD. Fig- 
ure shows that the widths of Mz>3 distributions for 
Xe+Sn at bombarding energies 32 to 65A MeV increase 
with the mean multiplicity according to a A = ^ scaling 
law, and this is confirmed by Figure Eb- The width for 
the 25A MeV system falls below this "systematic" trend, 
which may indicate that this energy is close to the thresh- 
old for multifragmentation in central collisions (the mean 
multiplicity for this system is approximately 3 fragments 
with Z > 3), leading to reduced fluctuations of the frag- 
ment multipHcity. Recent data obtained with INDRA for 
the same system at bombarding energies from 8 to 20A 
MeV will allow to study this point in more detail. 

Nevertheless the multiplicity of fragments does show 
some slightly more interesting features than Mtot- let us 
remark the "back-bending" in Figure for bombarding 



energies > 65A MeV, indicating the decrease of mean 
fragment multiplicity at the highest bombarding energies. 
The same behaviour is seen more clearly in Figure 03 for 
Ni+Ni collisions above 52A MeV. Let us note in passing 
that the available energy for the maximum of fragment 
production in the Ni+Ni system is here much lower (ISA 
MeV) than that reported in (17.5A MeV). 
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Figure 7: (Color online) Data for well-measured central colli- 
sions of Ni + Ni (6 < O.lbmax, Zdet > 0.8{Zproj + Ztarg)). (sl) 

Distributions of Mtot expressed in the variables of the second 
scaling law, i.e. Equation Q with A = i. The dashed curve 
refers to a global Gaussian fit to the data (dashed curve), (b) 
Log-log plot of the variance versus the squared mean value 
of the distribution of fragment multiplicity, Mz>3, for each 
bombarding energy. 

This type of "rise and fall" behaviour has been observed 
in many different data sets H, |^, Is^l- The absolute 
value of the maximum mean multiplicity and the en- 
ergy at which it occurs in central collisions are not only 
system-dependent |^ but also selection-dependent 
and detector-dependent. It should be noted that in the 
framework of the universal fluctuations theory this rise 
and fall is not consistent with the fragment multiplicity 
being an order parameter. Rather, in the energy domain 
for which fluctuations of Mz>3 show a regular scaHng 
behaviour (32-65 A MeV for Xe+Sn), it can be seen that 
this scaling is always of the second kind (A = i) and in 
this domain the multiplicity distributions all collapse to 
a unique, quasi-Gaussian scaHng function (Figure Eb). 

Therefore, neither the total charged multiplicity nor 
the fragment multipHcity have fluctuations which indi- 
cate the presence of different "phases" for central colli- 
sions in this energy range. It should be noted, however, 
that we cannot exclude the possibility of a different be- 
haviour of observables to which we do not have access 
in this data, for example the true total multipHcity in- 
cluding neutrons. Neither does our result exclude the 
possibility that in other fragmentation reactions, e.g. at 
higher bombarding energies and/or in spectator decays 
rather than central colHsions, the total or fragment mul- 
tiplicity may be the pertinent order parameter. For our 
data on central collisions between 25 and lOOA MeV how- 
ever we win from now on only consider the other possible 
order parameter for fragmentation, the size of the largest 
fragment. 
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Figure 8: (Color online) As Figure EK well- measured central 
Xe+Sn collisions) but for the charge Zmax of the heaviest 
fragment detected event-by-event, and A = 1. 



Looking at the log-log plot of the first two cumulant 
moments of Z^ax (Figure |HK) we can suspect some evo- 
lution of the scaling behaviour of this observable's fluc- 
tuations with increasing beam energy: the data do not 
appear to fall on a single straight line, but rather seem 
to be grouped into two "branches" with different slopes. 
Although most of the data in Figure |Ht he close to a Hue 
of slope A = 1, this does not correspond to a univer- 
sal scaling law, as all data in Figure IHb do not collapse 
onto a single universal curve under the flrst scaling law 
(Equation n] with A = 1). This confirms that there is a 
change in the scaling behaviour of Z^ax fiuctuations with 
increasing energy. The charge of the largest fragment 
in each event, and not the fragment multiplicity, seems 
therefore to be a good candidate for the order param- 
eter of multifragmentation in central collisions at these 
energies. 



C. Evolution of the scaling behaviour of order 
parameter fluctuations with incident energy 
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Figure 9: (Color online) (a) Zmax distributions for Xe+Sn 
collisions at 25-39A MeV, scaled according to Equation^with 
A given by a linear fit to the data of Figure |2K of the form 
a ~< m >^. The dashed curve is a best fit to scaled data 
using a Gaussian distribution, (b) As (a) but for bombarding 
energies 39-lOOA MeV. The dashed curve is a best fit to scaled 
data using the Gumbel distribution, Eq[Sl 



FigurelHt shows that the scaled Zmax distributions for 
Xe+Sn collisions from 45A MeV upwards are nearly iden- 
tical even in the large-Z^ax tails which are 2 orders of 
magnitude less probable than the most probable value of 
Zmax- The 39 A MeV data can also be included in this 
group if the small differences in the tail of the distribu- 
tion are neglected. However, the 32A MeV distribution is 
clearly significantly narrower, while the distribution for 
25A MeV is evidently of a different form, as can be seen 
both in the tails and around the maximum (see com- 
ments on the comparison of experimental distributions 
with this technique in Sec Hljl . 

Using the values of A extracted from FiglSK using a 
linear fit to the two 'branches' (25-39^ MeV and 39- 
100^ MeV) we find a good scaHng of distributions using 
a ~< Zmax >^l'^ at the lowest beam energies (FigureEt), 
whereas for the higher energies the fluctuations of the 
size of the largest fragment increase Hke cr ~< Zmax > 
(Figure Eb) ■ An approximately equally good scaHng for 
39A MeV data is achieved in both cases, and this energy 
may be considered as close to a transition between the 
two regimes. Figure El also shows that the shape of the 
scaling function changes with increasing beam energy: it 
is nearly Gaussian at low energy (the dashed curved in 
FigureEt represents a best fit to all data with a Gaussian 
distribution) , but at higher energies it is rather asymmet- 
ric with a near-exponential tail for Zmax greater than its 
most probable value (see Sec. IIV E|l . 

These observations establish the size (or charge) of the 
largest fragment as the most Hkely "order parameter" for 
fragment production in central collisions in this energy 
range. They also show that data can be assigned to one 
of two regimes depending on the A-scaling properties 
of this order parameter. At low energies systems obey 
the second-scaling law associated with an ordered phase, 
while at higher energies the first-scaling law is observed, 
typical of a disordered phase (large fiuctuations) It 
should be noted that these regimes are defined, not by 
the average size of the largest fragment produced for a 
given incident energy (which would be a typical way to 
define e.g. Hquid and gas phases), but rather by the way 
in which the fluctuations of this quantity evolve com- 
pared to its mean value. The identiflcation of the largest 
fragment size as order parameter also indicates some sim- 
ilarity between fragment production in central collisions 
and the aggregation models with the second-order phase 
transition mentioned in Sec. ^ whereas the monotonous 
behaviour of total and fragment multipHcities excludes 
"shattering" fragmentation models such as FIB (4^ . 

Another, not necessarily incompatible, interpretation 
of the different A regimes is suggested by Fig. ^ As we 
remarked in Sec. IIIIBI this figure suggests that at beam 
energies <39A MeV the heaviest fragment in central col- 
lisions is an evaporation residue of an (incomplete) fusion 
reaction. The data presented in Fig. for incident ener- 
gies 25-50A MeV correspond to the points highhghted by 
a circle in Fig. 31 Therefore the transition from a A ~ i 
to a A 1 branch can be interpreted as being linked 
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Figure 10: (Color online) As Figure |S| but for collisions of 
Ni+Ni from 32 to QOA MeV bombarding energy. 

to the disappearance of fusion-evaporation residues, and 
signalling either the onset of fusion-multifragmentation 
(phase transition) or the onset of transparency in central 
collisions (incomplete stopping). 

It should be noted that the transition energy of around 
39A MeV for the Xe+Sn system is slightly higher than 
that found in 01 (32^ MeV). This is due to the different 
method of selection of 'central coUisions'. In this paper 
we define central collisions based on the amount of en- 
ergy deposited into the light charged particles' degrees 
of freedom (i?ti2)- As we are dealing with finite sys- 
tems this reduces the energy available for the fragment 
degrees of freedom, due to energy conservation. This is 
consistent with the fact that when central collisions are 
selected based on fragment degrees of freedom (the frag- 
ment kinetic energy fiow angle selection of 0, 0, Ef) 
the multifragmentation regime is observed at lower inci- 
dent energy for the same system (32^ MeV for Xe+Sn 
in H^), as is the transition to the A = 1 scaHng regime. 

The observed scaling properties of Z^ax fiuctuations 
are confirmed by the Ni+Ni data, as shown by Figure EH 
and Figure El Most data follow a A « 1 scaHng law 
for Zmax fiuctuations, except at the lowest energies. The 
data do not all collapse to a single scaling function in 
terms of a A = 1 law (Figure llOb l , but only for bom- 
barding energies E > 52 A MeV (Figure Eb)- The simi- 
larity between the scaling functions observed for the two 
different systems should be noted, both in the "ordered" 
(A - i) and "disordered" (A - 1) regimes (Figures 
inland Table mi. 



D. System-size dependence of energy of transition 
from "ordered" to "disordered" regime 

We observe for the Xe+Sn and Ni+Ni data that the 
energy ranges corresponding to the different regimes are 
not the same for two systems of different total mass, the 
"transition" occurring around 39A MeV for Xe+Sn and 
52A MeV for Ni+Ni. The "disordered" regime begins at 
lower incident (or available) energy for the heavier sys- 
tem. This tendency is confirmed by the data for Ar+KCl 
and Au+ Au (figures El and I13|l . The former has a total 



Figure 11: (Color online) (a) Zmax distributions for Ni+Ni 
collisions at S2-52A MeV, scaled according to EquationQwith 
A given by a linear fit of the form a ~< m >'^. The dashed 
line is a best fit using a Gaussian distribution, (b) As (a) but 
for energies 52-90A MeV. The dashed line is a best fit using 
the Gumbel distribution, Eq[^ 
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Figure 12: (Color online) As Figure El but for collisions of 
Ar+KCl from 32 to TiA MeV bombarding energy. 



mass smaller than both Xe+Sn and Ni+Ni. The largest 
fragments produced in central coUisions of ■^^ Ar+KCl fol- 
low quite closely a A '--^ i scaling behaviour for beam 
energies up to at least 52 A MeV (Figure 112b.). while 
the data for 74A MeV, which deviate from the overall 
trend, may indicate a transition to the A ~ 1 regime 
occurs somewhere between these last two available data 
points. The data for 32-52A MeV exhibit very similar 
near-Gaussian probability distributions. Although we 
can not have much confidence in a scaling law estab- 
lished for only two data points let us remark in passing 
that the scaling behaviour observed for data at 52 and 
74A MeV is approximately a A 1 scaHng with a scaHng 
function of a form similar to that observed for the other 
"disordered" regime data, indicating that the transition 
energy for this system is probably somewhere between 52 
and 74^1 MeV, higher than for the heavier systems. 

On the other hand, for the much heavier ^^^Au+^^^Au 
system a A = 1 scaHng law is observed for fiuctuations 
of the size of the largest fragment in each event for all 
studied beam energies (Figure I13|l . Even at the lowest 
energy (40A MeV) the scaHng of the entire Z^ax dis- 
tribution with the other energies is exceHent. For this 
system it may be that the beam energy of the transition 
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from the A = i regime, if it exists, to A = 1 occurs at a 
lower energy than the lowest energy available for study 
in our data, which may explain why we do not observe an 
"ordered" regime for this system. Of course it is equally 
possible that no such regime exists for this system. In- 
deed, following the interpretation of Sec lIV CI in terms 
of the disappearance of fusion-evaporation residues, we 
would not expect fusion to occur for a system as heavy 
as Au+Au at any incident energy. 
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Figure 13: (Color online) As Figure |H1 but for collisions of 
Au+Au from 40 to 80^4 MeV bombarding energy. 

We therefore observe that the bombarding energy (or 
available energy) at which there is a transition from the 
"ordered" to the "disordered" regime decreases with in- 
creasing total mass of the system under study. This is 
represented in Figure If one interprets our results in 
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Figure 14: Dependence on bombarding energy and total sys- 
tem mass of the frontier between the two A-scaling regimes 
observed in this work for very central collisions. 

the framework of the universal fluctuations theory, i.e. 
in terms of a 2nd order phase transition, one would ex- 
pect the energy at which one phase is replaced by the 
other to be related to the critical temperature of the cor- 
responding system. The deflnition of this quantity for 
flnite, charged, systems such as nuclei has received much 
theoretical attention. On the one hand, in the absence 
of Coulomb forces, a pseudo-critical temperature (cor- 
responding to a large peak in the flnite system speciflc 



heat) may be deflned which is smaller than the inflnite 
matter value Tc , and increases as the size of the system 
approaches the inflnite matter limit [fi^l- On the other 
hand, calculations including the Coulomb repulsion show 
that the maximum temperature that an equilibrated hot 
nucleus can support , Tum, decreases for heavier nuclei, 
due to their increasing nuclear charge [63] . A recent sys- 
tematic study of a wide range of data on so-called caloric 
curves by Natowitz et al |^ has shown that the tem- 
perature and excitation energy at which a plateau is ob- 
served in these curves decreases with increasing mass of 
the primary excited nucleus, and that such behaviour 
is consistent with theoretical predictions for the limiting 
temperatures, Tum- 

In the present analysis we chose not to study this ques- 
tion in more detail, as to do so would require: (i) to 
show that the fragment production is thermally-driven; 
(ii) to identify the thermal (sub-)system in each case; 
and (iii) to deduce the mass, charge, excitation energy 
and temperature of this (sub-)system. Each step would 
require us to make important hypotheses about the for- 
mation and decay properties of fragments in the colli- 
sions under study, taking us far from our initial goal of 
determining as much information as possible on the na- 
ture of fragment production with a minimum number 
of suppositions. Moreover, the observed effect, which 
depends on the entrance channel total mass and avail- 
able/bombarding energy may have a completely different 
origin. Therefore we will limit ourselves to the observa- 
tion that the transition from "ordered" to "disordered" 
regime takes place at a lower available energy for sys- 
tems of greater total mass. 



E. Detailed study of the form of the scaling 
functions 



The results of the analysis for all the systems pre- 
sented here are summarized in Table Systems have 
been grouped according to their observed A-scaling be- 
haviour, A ~ i at the lowest beam energies and A 1 
at the highest. In order to get some quantitative infor- 
mation on the form of the scaling function in each case, 
we calculated the overall coefflcient of skewness, 7, and 
the kurtosis, k, for each system in each regime, using the 
following deflnitions 
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Larger skewness values indicate more asymmetric distri- 
butions, while the kurtosis measures the deviation of the 
distribution from the Gaussian form (k = being a per- 
fect Gaussian). The values shown in the table conflrm 
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Table I: Fluctuation scaling exponent A, coeflEicient of skewness 7, kurtosis k, ^ for global fits to data with Gaussian (A ~ ^) 
and Gumbel (A ~ 1) distributions, and scaling function tail exponent f'(x^) (see Section^, for the scaling functions shown in 
Figures l9l lTT1IT2l andlT3l 



System 


£;,„c (AMeV) 


A 


7 








Xe+Sn 


25-39 


0.48 


0.45 


0.66 


27 




Ni+Ni 


32-52 


0.61 


0.48 


0.64 


72 




Ar+KCl 


32-52+ 


0.46 


0.33 


0.64 


63 




Au+Au 


40-80 


1.00 


0.83 


1.39 


1.8 


0.92±.09(1.3) 


Xe+Sn 


39-100 


0.89 


0.85 


1.20 


11 


1.20±.03(1.2) 


Ni+Ni 


52-90 


0.88 


0.84 


1.30 


45 


1.54±.04(7) 



our observation that the higher-energy, A ^ 1-scaUng 
data present more asymmetric, less Gaussian probabil- 
ity distributions than those at lower energy. It can also 
be seen that, quantitatively, the scaling functions cor- 
responding to the two regimes for systems of different 
masses have very similar forms, confirming the fact that 
data for different systems do indeed collapse on to a sin- 
gle distribution. 

In the "ordered" regime the fluctuations of the size of 
the largest fragment show a significant deviation from a 
Gaussian distribution, on the contrary to the near-perfect 
Gaussian distributions that we observe for total and frag- 
ment multiplicities (Figures El El and . For the "disor- 
dered" regime the deviation from the Gaussian form is 
large and the shape of the distribution function is well 
reproduced by a Gumbel distribution, 

$(za) ~ exp- (z - e"^) (5) 

This is shown by the dashed curves in Figures IHb and 
ITTh . and the grey curve in Figure [T3h. which represent 
the best fits to these data using Equation [3 The overall 
agreement can be seen to be quite good, and excellent for 
the Au+Au data where = 1-8 is achieved. Some sig- 
nificant deviations can be seen however in the tail of the 
distributions for the lighter systems Ni+Ni and Xe+Sn: 
the data seem to have a faster-than-exponential fall-off 
for large (positive) fiuctuations about the mean value. 

To confirm this in a quantitative manner we performed 
fits to the tails of these functions with the asymptotic 
scaling function form exp— z'' (see Table Pi . They show 
that the exponent v is very close to 1 for the Au+Au data, 
confirming the observation of an exponential tail compat- 
ible with the Gumbel distribution, while for Ni+Ni and 
Xe+Sn data we find 1 < D < 2. Let us remark in passing 
that, just as A seems to increase toward the asymptotic 
value of 1 with increasing system mass in the "disordered" 
regime, it is possible that v also has a systematic mass 
dependence, and decreases towards an asymptotic value 
of 1 (exponential tail) for the heaviest system. 

We have clearly and quantitatively established the 
form of the scaling functions in the two regimes, and in 
the data presented here we do not observe any deviation 
from these "canonical" forms, quasi-Gaussian at low en- 
ergy and quasi-Gumbel at high energy. We do not, there- 



fore, have any information on the nature of the transition 
between the two regimes: as we discussed in Section UTI 
we cannot be certain of the form of the scaling function 
at the critical point, but one may see a sharp decrease 
of the large-z tails (faster than exp — z^), or an order pa- 
rameter distribution (OPD) with an exponential large-z 
tail but significantly different from the Gumbel distribu- 
tion around the maximum and below On the other 
hand, for a first-order phase transition with a passage 
through the coexistence region, although the sca ling be- 
haviour of the OPD is not well established ^2, by 
definition in this case the order parameter should present 
a bimodal (double-humped or very wide) distribution. 
More detailed study of existing and new data around the 
transition energy of 39A MeV (for the Xe+Sn system) 
may yet reveal such features. 



F. Dependence of the observed scaling behaviour 
on the violence/centrality of collisions 




ln(K,2=<m>^) z,4, 

Figure 15: (Color online) As Figure |2l but for collisions with 
an estimated centrality of & < 0.2bmax- 

The results presented above come from a wide-ranging 
set of data concerning different system energies and 
masses, but they represent only 1% of the total mea- 
sured cross-section due to our centrality cut. The repre- 
sentativity of the observed scaling behaviour is far from 
certain in this case, as it is well known that data se- 
lected by cuts in distribution tails are prone to serious 
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autocorrelation effects due to conservation laws (iH l63|. 
We therefore have to study the dependence of our find- 
ings on the strictness of our centrality cuts. In Figure [151 
the data for the same Xe+Sn collisions as in Figure 
have been analyzed with a slightly relaxed centrality cut, 
b < 0.2bjnax- The effect on the apparent value of A can 
be seen from the cumulant moments plot, Figure lT^h ,: re- 
laxing the centrality condition gives an apparent scaling 
law with a smaller value of A, which is here A « 0.85 in- 
stead of A = 0.89 (see Tabled. A similar dependence on 
the estimated centrality of collisions was observed in 0|, 
and interpreted in terms of a smaller excitation energy 
of fragmenting systems in less central collisions. How- 
ever, as Figure [T5h shows, the A-scaling actually begins 
to break down for these less central data, as can be seen 
in the high-Z(A) tail of the "scaling" function. 
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Figure 16: (Color online) As FignrelT^ but for collisions with 
an estimated centrality of 6 < ^Ahmax- 

If the centrality condition is relaxed even further, 
h < OAb„iax (Figure El, we first observe that now ap- 
parently all of the data for Xe+Sn is compatible with 
an approximate second-scaling law (A = 0.58), without 
any change of scaling regime between 25 and lOOA MeV. 
However we can not assign the data to a single "ordered" 
regime because in fact this is not true A-scaling: the 
different distributions no longer collapse to a universal 
curve (Figure Eb)- 

Therefore we observe universal fluctuations and be- 
haviour of the Zmax observable compatible with it be- 
ing an order parameter for nuclear multifragmentation 
only for very central collisions of symmetric systems for 
which one may suppose a near-to-total overlap of the 
projectile and target in the entrance channel. We have 
not, up to now, observed an equivalent scaling for quasi- 
projectile residues in mid-peripheral to peripheral colli- 
sions. This may be because in such reactions the frag- 
ment production is far more sensitive to entrance channel 
effects, which vary greatly as a function of the colliding 
nuclei and their energies. If so, a meaningful comparison 
between different systems is harder to achieve, at least 
with our rather "global" approach. In head-on collisions, 
on the other hand, the fragmentation of the projectile 
and target may be virtually a statistical (although not 
necessarily thermal) process, far less sensitive to the de- 
tails of the reaction, and therefore more amenable to re- 



veal features which are independent of the system under 
study. 



V. SUMMARY 

We have studied nuclear multifragmentation data ob- 
tained with the INDRA 4tt array for collisions of symmet- 
ric systems of total mass Atot ~ 75 — 400 at bombarding 
energies from 25 to lOOA MeV. Using the total transverse 
energy of light charged particles, £'ti2, as a measure of 
collision violence, we deduced the evolution of these re- 
actions with beam energy, impact parameter and system 
size from the experimentally measured correlations be- 
tween the charge of the largest fragment detected in each 
event, Zmax, and £"412. For all data presented in this 
work these correlations are dominated by reactions lead- 
ing to a projectile-like fragment whose size decreases with 
increasing collision violence. In the case of the heavi- 
est projectiles (Xe, Au) fission of the moderately excited 
quasi-projectile modify this picture. For higher excita- 
tion energies the opening of the quasi-projectile multi- 
fragmentation channel may exhibit a bimodal behaviour 
which has been evidenced [g^ in selecting complete de- 
tection events. On the other hand, in central collisions 
for a few percent of the measured cross-section "heavy" 
residues are produced in the systems Ar+KCl, Ni+Ni 
and Xe+Sn, suggesting incomplete fusion of projectile 
and target, for beam energies which are not too high 
(< 52A MeV for Ar+KCl, < 40^1 MeV for Ni+Ni, < 32A 
MeV for Xe+Sn). 

Using a model-independent analysis based on the the- 
ory of universal fluctuations of the order parameter for 
flnite systems, we tested the most violent collisions for 
signals that the fragment production may be related to 
a phase transition. Following the results of we first 
confirmed that, of the two possible order parameters for 
a critical fragmentation process, it is the charge (size) 
of the largest fragment, Zmax, and not the total or IMF 
multiplicities Mtot or Mz>3, which has a behaviour of 
the scaling properties of its fiuctuations compatible with 
its being an order parameter for a critical fragmenta- 
tion process. Indeed, we have shown that the event-by- 
event distribution of Zmax allows to sort data into two 
fiuctuation-scaling regimes defined by the value of the 
scaling exponent A, which is approximately equal to ^ 
at low energies and tends towards the asymptotic value 1 
at high energies with increasing total system size. These 
regimes are equally well characterised by a distinctive 
form of the scaling function ^{z(/^)). At low energies 
this function, although more symmetric than that seen 
at higher energy, is significantly different from the Gaus- 
sian form. The deviation is quantitatively the same for 
the three different-sized systems (Ar+KCl, Ni+Ni and 
Xe+Sn) for which we observed this low-energy regime. 
In the high-energy regime the scaling function is more 
asymmetric and tends towards the asymptotic form of the 
Gumbel distribution with increasing system mass. This 
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evolution concerns mainly the \arge-Zmax tail of the dis- 
tribution which falls off more slowly for heavier systems, 
becoming exponential. For the Au+Au system $(2(a)) 
is an almost perfect Gumbel distribution. 

The bombarding energy at which the passage from one 
regime to the other is situated decreases as the total sys- 
tem mass (and charge) increases. This is the contrary 
of the expected behaviour if this transition were related 
to the critical temperature of the systems under study, 
which would increase with the mass, but is on the other 
hand consistent with the observation of decreasing limit- 
ing temperatures for finite nuclei |3| due to the increase 
in Coulomb energy for heavy nuclei. However, we can 
only speculate whether this is the reason for our obser- 
vation, as our analysis does not depend on any model 
of fragment production such as supposing it to be equili- 
brated or thermally-driven, and therefore it does not give 
any information on whether these conditions are met or 
not. On the other hand, data strongly suggest that the 
dependence on entrance channel mass and bombarding 
energy of the two regimes is closely linked to the disap- 
pearance of heavy residues in central collisions. 

Although this analysis allows to establish the existence 
of two distinct regimes in multifragmentation reactions 
based on the scaling properties of the fluctuations of the 
largest fragment, it has not given any information on 
the passage from one regime to the other. In the data 
studied so far we observe only order parameter distribu- 
tions which are compatible with weakly-correlated sys- 



tems, i.e. far from the critical point. It may be that the 
general survey of a wide range of data presented here is 
not sufficiently detailed to reveal such features, or that 
data taken in small bombarding energy steps around the 
"transition energy" is necessary to track the evolution of 
the order parameter distribution. Such additional data 
for the Xe+Sn system have been measured recently and 
are currently under analysis. It is our feeling that the 
present work constitutes a solid basis for further study of 
the question of criticality in nuclear multifragmentation 
data. 
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